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'  An  analytic  study  is  made  of  the  effects  produced  by  a  localized 
modulation  of  the  electrical  conductivity  of  a  magnetized  plasma  carrying  a 
field  aligned  current.  The  modulation  is  assumed  to  be  caused  by  a  heating 

source  external  to  the  plasma.  The  formulation  used  is  applicable  to 

•  >  •  ,  U-  •  "S’  O  Id  t  v.  » 

modulation  frequencies below  the  ion  cyclotron  f requency^to^ ^  and  permits 
the  calculation  of  the  radiation  pattern  which  contains  shear  as  well  as 
compressional  Alfven  modes,  and  near  fields.  For  u/u>cj^<<  1,  the  radiation 
appears  primarily  in  the  form  of  a  highly  collimated  shear  mode  travelling 
along  the  field  lines  passing  through  the  spot  where  the  modulation  occurs. 
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l.  Introduction 

A  possible  method  of  exciting  low  frequency  electromagnetic  waves  consists 
of  the  external  remote  modulation  of  zero  order  current  systems  flowing  In 
magnetized  plasmas.  The  motivation  for  taping  internal  currents  Is  that  the 
dimensions  of  external  antennas  required  for  low  frequencies  are  prohibitively 
large.  A  significant  shortening  of  the  effective  length  of  antennas  formed 
inside  a  plasma  results  from  the  density  dependence  of  the  Alfven  speed,  i.e., 
VA  =  B0/(4irMN) 1 /2  where  B0  is  the  confining  magnetic  field,  M  is  the  ion  mass, 
and  N  the  plasma  density.  The  characteristic  dimensions  of  an  antenna  that 
excites  an  Alfven  wave  in-situ  is  shortened  from  the  vacuum  value  by  a  factor 
nA  *  C^VA  “  a)pi/(*)ci»  where  n^  is  the  Alfven  index  of  refraction,  u)p^  is  the 
ion  plasma  frequency,  is  the  ion  cyclotron  frequency,  and  c  is  the  speed 
of  light.  Over  a  broad  class  of  plasmas,  both  in  the  laboratory  and  naturally 
occurring,  this  shrinkage  ranges  from  a  factor  of  10  to  103. 

Although  there  is  a  clear  attractiveness  in  this  method  with  regard  to 
antenna  size,  the  radiation  efficiency  and  directionality  of  the  radiated 
signals  must  be  carefully  examined  in  order  to  assess  the  overall  usefulness 
of  this  scheme.  It  is  the  purpose  of  the  present  analytical  study  to  address 
these  issues  by  calculating  the  radiation  pattern  resulting  from  the  thermal 
modulation  of  a  current  system  frequently  found  in  plasmas. 

Specifically,  this  study  considers  the  external  modulation  of  the  elec¬ 
trical  conductivity  of  a  collisional  plasma  carrying  a  steady  ohmic  current 
along  the  confining  magnetic  field,  i.e.,  a  field  aligned  current.  The 
thermal  modulation  is  assumed  to  be  periodic  with  frequeny  co  <  u>c:j ;  It  Is 
envisioned  that  the  thermal  modulation  arises  due  to  a  high  frequency  electron 
heating  source  with  frequency  <*>rf  >>  vg>  where  vg  la  the  effective  electron 


collision  frequency 


In  practice,  u)rf  would  be  near  one  of  the  typical  electron  RF  heating  frequen¬ 
cies,  i.e.,  the  electron  plasma  frequency  Wpe  or  electron  cyclotron  frequency 
wce.  In  the  present  study  we  do  not  address  the  issue  of  thermal  transport 
associated  with  the  heater  wave  at  u)rf.  It  is  simply  assumed  that  a  localized 
hot  spot  of  Gaussian  shape  is  formed  at  some  point  in  an  otherwise  uniform 
plasma.  This  choice  of  a  profile  simplifies  the  analysis  and  is  qualitatively 
consistent  with  a  detailed  heat  transport  study  previously  performed  by  these 
authors  for  the  high  latitude  ionosphere.1 

One  of  the  features  of  the  present  analysis  is  that  in  calculating  the 
radiation  of  Alfven  waves  no  expansion  in  t»/u>cj  is  made.  This  implies  that 
the  relative  efficiency  for  radiating  the  shear  Alfven  mode  and  the  compres- 
sional  Alfven  mode  (which  are  simultaneoulsy  generated)  is  found.  For 

w/u >ci  “0.2  nearly  90%  of  the  radiated  power  appears  in  the  form  of  a  shear 

Alfven  mode  which  travels  along  the  magnetic  flux  tube  passing  through  the  hot 

spot.  The  compressional  Alfven  mode  carries  away  the  remaining  10%  in  a 

quadrupole-like  pattern  that  moves  out  of  the  flux  tube  containing  the  shear 
mode.  The  efficiency  of  the  process  is  of  the  order  of  10-5,  which  is  quite 
low  on  absolute  standards,  but  may  result  in  detectable  signals  for  the 
auroral  ionosphere. 

The  principal  results  of  the  study  are  as  follows.  The  radiation  pattern 
due  to  modulation  of  the  electrical  conductivity  is  found  for  u)/wcj  <  1.  It 
contains  a  shear  wave,  a  compressional  wave,  and  near  fields.  The  shear 
Alfven  mode  dominates  and  results  in  a  headlight  radiation  pattern.  The  low 
frequency  highly  directive  wave  electric  fields  are  of  the  order  of  (ATe/Te)EQ 
where  E0  is  the  zero  order  dc  field  driving  the  field  aligned  current,  and 
ATe/Te  is  the  fractional  modulation  of  the  electron  temperature. 

The  manuscript  is  organized  as  follows.  In  Sec.  II,  the  plasma  model  is 


described  and  the  radiation  fields  are  calculated.  Various  parameter  depen¬ 
dences  of  the  radiation  pattern  are  discussed  in  Sec.  Ill,  and  conclusions  are 
presented  In  Sec.  IV. 

II.  Analysis 

Consider  an  Infinite,  uniform,  colllsional  plasma  In  which  a  zero  order 

A  A 

steady  current  density  J0z  flows  along  a  uniform  magnetic  field  B0z.  An 
external  high  frequency  source  of  carrier  frequency  u>rf  is  amplitude  modulated 
at  a  low  frequency  to  such  that  the  following  frequency  ordering  holds 

vi  «  “  <  “cl  «  ve  «  “pi  «  “rf  ~  "ce.  “pe  <*> 

In  here  Vj  represents  the  ion  collision  frequency,  and  v0  *  vei(Te)  +  ven^Te^ 
is  the  total  effective  collision  frequency  for  electrons  due  to  collisions 
with  ions  (vel)  and  with  neutrals  (ven),  both  of  which  depend  on  the  electron 
temperature  Te.  The  amplitude  modulated  high  frequency  source  is  assumed  to 
produce  a  periodic  variation  in  Te  whose  magnitude  enters  into  the  present 
study  as  an  input,  but  which  in  general  must  be  obtained  from  an  appropriate 
heat  transport  calculation.  The  Te  modulation  is  assumed  to  be  spatially 
localized  with  a  Gaussian  shape  having  widths  d|  and  dj^  In  the  directions 
parallel  and  perpendicular  to  the  magnetic  field,  respectively.  Spatial 
localization  of  the  modulation  In  Te  can  in  practice  be  brought  about  by 
focusing  the  high  frequency  source  through  appropriate  external  antennas  or  by 
using  natural  resonances  or  cut  offs  of  the  plasma  (e.g,  u)rf  ~  Wpe,  u>ce). 

The  sinusoidal  temporal  variation  In  Te  gives  rise  to  an  equivalent 
variation  in  the  electrical  conductivity  tensor.  For  the  ordering  indicated 
in  Eq.  (1),  the  principal  effect  is  to  change  the  component  parallel  to  the 
magnetic  field  giving  Aoz  ~  -ozo(Ave/ve)  ~  ozo(ATe/Te).  We  express  this 
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where  e  and  m  refer  to  the  electron  charge  and  mass,  E0  Is  the  zero  order  dc 

electric  field  responsible  for  sustaining  the  field  aligned  current,  and 

o(u>)  =  a  e  +  Oj  is  the  complex  conductivity  tensor  due  to  electron  (oe)  and 
«• 

ton  (0j)  contributions.  In  Eq.(3)  refers  to  the  self-consistent  electric 
field  at  frequency  generated  by  the  plasma  In  response  to  the  change  in 
temperature  ATe.  The  effect  of  appears  in  two  forms:  1)  through  the  usual 

dielectric  response  represented  by  o(u>)  =  (iu)/4it)(I  -  e)  with  e  the  dielectric 
tensor;  2)  as  a  space  charge  bunching  due  to  the  pile  up  of  flowing  electrons, 
and  represented  by  the  term.  The  third  term  in  Eq.(3)  is  the  effect  of 
the  thermo-electric  modulation. 

We  shall  show  that  the  last  two  terms  In  Eq.(3)  can  be  made  negligible 


ion  our  problem  and  we  will  only  use  the  current  due  to  the  ohmic  modulation 
as  source.  Eq.(3)  is  then  rewritten  as  the  sum  of  the  wave  dielectric  and  the 
source  contributions 

J.  =  o(to)*_E  +  Aoz  E0z  .  (3a) 

Since  the  charge  bunching  term  arises  due  to  motion  parallel  to  the  mag¬ 
netic  field,  its  effect  is  identical  to  that  in  an  unmagnetized  plasma,  i.e., 
it  contributes  a  correction  to  the  Aoz eoz  source  in  Eq.  (3)  of  order  (l6/lt)2, 
where  *  (c2/4iru)Ozo)  1  /2  is  the  collisional  diffusion  length  (or  skin  depth) 
and  is  the  spatial  scale  length  over  which  the  temperature  modulation  occurs. 
To  obtain  coupling  to  an  Alfven  wave  one  expects  that  £t  ~  2x/kA  =  2ttva/u>; 
this  implies  that  the  space  charge  term  contributes  a  correction  of  order 
(m/M)(u>ve/iocj2).  In  the  present  study,  parameters  are  chosen  so  that  this 
correction  plays  a  negligible  role.  The  thermo-electric  modulation  arising 
from  the  ATe  terra  in  Eq.  (3)  is  present  even  when  E0  *  0,  i.e.,  in  the  absence 
of  field  aligned  currents.  This  term  can  also  contribute  to  the  radiation  of 
Alfven  waves  due  to  both  the  parallel  and  perpendicular  components  of  the 
gradient.  The  parallel  contribution  becomes  relatively  small  for  E0»Te/e£t, 
while  the  perpendicular  contribution  ceases  to  dominate  for  EQ»  ( u>cj/w)Te/e£t. 
The  amplitude  of  E0  is  limited,  however,  because  in  order  for  the  steady-state 
ohmic  description  to  be  applicable,  the  value  of  E0  must  be  smaller  than  the 
Dreicer  runaway  field.  These  constraints  imply  that  in  order  for  the  field 
aligned  currents  to  play  a  significant  role  in  the  radiation  of  Alfven  waves 
the  following  inequalities 

V1  »  <eEQ/Te)  »  J^-1  (4) 

must  be  satisfied.  In  Eq.  (4)  lc  =  v/vg  is  the  collisional  mean  free  path  for 
electrons,  which  of  course,  must  satisfy  lc  «  lt  ~  vA/u>  and  v  is  the  electron 


thermal  velocity.  In  the  present  study  we  assume  that  the  system  under  consi 
deration  satisfies  Eq.  (4)  since  otherwise  the  effect  of  field  aligned  cur¬ 
rents  would  be  masked  by  modulations  accompanying  the  thermo-electric  fields. 

Next  we  proceed  to  determine  the  electric  field  _E.  We  Fourier 
analyse  in  configuration  space  and  obtain 


E(k)  = 


exp[-i(kxx  +  kyy  +  kzz)]  E(£> 


(5) 


Maxwell's  wave  equation  yields 

k  x  k  x  E  =  -  k02(E  -  -7^-  J)  ,  (6) 

—  —  —  —  1 0)  — 


where  is  the  Fourier  transform  of  and  kQ  =  w/c.  From  F,q.(3a),  the 


Fourier  transformed  current  of  Eq.(6)  may  be  written  as 


J  *  o(o))*i  +  zozo  AF(kx,k  kz)  ,  A  =  (— — )  Eq  ,  (7) 

vzo/  m 


where  F  is  the  Fourier  transform  of  f(x,y,z),  and  where  we  have  assumed  that 
the  parallel  ac  conductivity  oz  -  ozo,  the  unperturbed  dc  value.  This 
assumption  is  correct  for  the  frequency  range  specified  in  Eq.(l). 

We  solve  Eq.(6)  for  the  unknown  field  j£  and  obtain 


E 


x,y 


kz  &x,y 
DC jc ,  u>) 


AF 


(«) 


with  the  notation  Ex>v  meaning  either  Ex  or  Ey  and  likewise  for  other 
variables.  In  Eq.  (8)  the  coupling  factor  is 

gx,y  =  kx,y  (k2  “  ko2  ^l)  ±  ky,x  ko2  exy 


(9) 


I 
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with  the  (+)  sign  for  gx  and  the  (-)  sign  for  gy.  In  here  is  the  perpen¬ 
dicular  component  of  the  dielectric  tensor  for  a  cold  collisional  plasma  and 
exy  the  corresponding  off-diagonal  component,  the  parallel  component  is  £|. 

The  denominator  in  Eq.  (8)  leads  to  the  dispersion  relation  containing  both 
the  shear  and  compressional  Alfven  modes  and  is  more  clearly  expressed  in  the 
form 

D(k,u>)  =  (kz2  -  ks2)  (kz2  -  kc2)  ,  (10) 

with 

k2s,c  -  C  eI  1  0/2)  Uki4  (1  -  ei/e,)2  +  4k04  £xy  cyx  (1  -  ki2/k02  e,)]1/2 

+  kX2  (1  +  EJ./G,)}  ,  (11) 

where  the  upper  sign  goes  with  ks 2  and  the  lower  sign  with  kc2.  Physically, 
the  kg  and  kc  roots  correspond  to  the  parallel  wave  number  of  the  shear  and 
compressional  modes,  respectively,  for  fixed  values  of  kj^  »  (kx2  +  ky2)1/2 
and  u>.  For  the  sake  of  clarity  the  characteristic  topology  of  these  roots  for 
a  fixed  value  of  u  <  a>cl  is  shown  in  Fig.  2  as  a  function  of  k^2.  As  is  well 
known,  the  shear  mode  is  essentially  insensitive  to  variations  in  k^,  while 
the  compressional  mode  exhibits  an  abrupt  cut-off  for  a  k^  slightly  below  kA. 
The  topology  of  these  roots  plays  an  important  role  in  the  evaluation  of  the 
radiation  pattern  associated  with  the  hot  spot.  It  should  be  noted  that  in 
the  limit  u>/u>ci  <<  1,  kg2  -  kA2  (1  +  iVj/u>)  and  kc2  =  kA2  (1  +  iVj/u>)  -  k^2. 


4 


fl 


Because  of  the  high  electron  conductivity  along  the  magnetic  field,  ej  is  a 
large  quantity  having  negligible  contribution  on  low  frequency  waves  as  may  be 
clearly  seen  in  F.qs.  (8)  and  (11)  when  we  let  £||  ■*■  However,  it  is  kept  in 
the  analysis  for  the  sake  of  completeness.  The  large  E|j  approximation  is  also 
used  in  calculating  the  z  component  of  the  electric  field,  which  is  obtained 
either  in  k  space  form  Eqs.  (6)  and  (7)  as 


Ez  =  -  AF(kx,  ky ,  kz) 


or  in  configuration  space  from  the  z  component  of  Eq.(3a) 


/Aoz  (x,y,z)\ 

=  v  ^r~r°  • 


(12a) 


where  the  contribution  from  Jz/aZQ  is  negligibly  small.  Current  modulation 
parallel  to  the  magnetic  field  is  therefore  negligible,  with  Ez  keeping  the 

A 

Initial  current  -I0z  nearly  constant.  The  transverse  variation  of  Ez  implies 
that  transverse  electric  fields  are  generated  and  associated  polarization 
currents  are  forced  to  flow  in  the  region  where  the  hot  spot  is  formed.  These 
locally  forced  currents  create  a  self-consistent  field  in  the  neighborhood  of 
the  hot  spot  which  appears  in  the  analysis  in  the  form  of  a  near  field.  In 
addition,  outwardly  propagating  Alfven  waves  emanate  from  the  locally  forced 
currents;  these  waves  have  quasi-transverse  electric  fields  with  Ez  negli¬ 
gibly  small  outside  the  hot  spot. 

Next  we  proceed  to  Invert  the  Fourier  transformed  field  _E  to  obtain  the 
radiation  pattern  in  configuration  space.  Since  the  source  has  azimuthal 
symmetry,  it  is  advantageous  to  use  a  cylindrical  coordinate  system  in  k  space 
as  well  as  in  configuration  space  to  obtain 


» 
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Ex,y(P,  ♦,  2)  =  — 


00  00 

^3  J* lkl  R  jd 


00  2tt 

C  ^zSx  v^l  ^ 
k,  Ida  F(k1,kz) 

2  /  (k  2-k  2)(k  2-k  2) 


exp(i[kzz  +  kj^P  cos(a-<j>)]} 


where  (p,<j>)  are  the  radial  and  azimuthal  coordinates  In  configuration  space 
and  (k^,a)  are  the  conjugate  quantities  In  k.  space. 

It  should  be  noted  that  In  the  Integrand  of  Eq.  (13)  the  only  a  depen¬ 
dence,  aside  from  the  expected  exponential  factor,  arises  through  gx>y  and 
that  such  a  dependence,  as  seen  from  Eq.  (9),  is  proportional  to  either  sina 
or  cosa.  This  implies  that  the  required  Integrations  over  a  are  of  the  form 


da(slna, cosa)exp [  1  k^p  cos(a-<|>)]  =  2ui  Ji(k^p)(sin  <|>,cos  j>),  (14) 


where  Jj  Is  the  first  order  Bessel  function. 

Expressing  the  electric  field  components  in  cylindrical  coordinates 
through  the  relation  Ex  =  Ep  cos<J>  -  E^  sim}>,  Ey  =  Ep  sin<f>  +  E^  cos4>  and  using 


Eq.  (14)  yields, 


En  <f)  (  p  )^|Z) 


-1A  [1  .  2  ( 1  rF(kl-k»>Jl 

(2  «>2J  1  1  J  Z  <kz2-k,2 


[1(k1p)Gikzz  L 

,2>(kz2-kc2)  “•*  •  "5> 


where  again,  the  factor  hp  -  (kj_  +  kz  -  k0  e^)  goes  with  Ep  and  h^  = 
(~k02  exy^  wlth  E<t>»  and 


E(ki,kz)  -  ( tt  ) 3  / 2  (d||  d_L2)  exp  {-  ^  [kx2  dj_2  +  kz2  dp2]}. 


(16) 


Next  we  proceed  to  perform  the  kz  integration  by  making  a  suitable  dis¬ 


tortion  of  the  contour  of  integration  in  the  complex  k7  plane.  For  z  >  0 
the  relevant  features  of  the  integrand  are  the  existence  of  two  simple  poles 
at  ks  and  kc  with  Im  (ks,  kcJ  >  0,  and  a  saddle  point  at  kSp  =  i2z/d||2,  as  is 
indicated  in  Fig.  3.  It  is  evident  from  Eq.  (11)  that  the  location  of  the 

poles  is  a  function  of  kx,  which  is  also  shown  in  Fig.  3.  For  kx  0  the  two 

poles  are  close  to  each  other,  and  as  kx  increases  the  kc  pole  goes  to  infinity 
along  a  line  parallel  to  the  imaginary  axis.  The  ks  pole,  however,  remains 
essentially  fixed  as  kx  varies.  The  saddle  point  is  independent  of  kx  and  its 
noteworthy  feature  is  that  for  some  value  of  kx  the  kc  pole  may  be  close  to 
it.  Such  a  behavior  can  be  handled  by  applying  a  technique  prescribed  by 
Banos  and  outlined  in  Appendix  A.  Since  the  hot  spot  is  symmetric  in  z, 

analogous  features  in  the  lower  part  of  the  complex  kz  plane  are  relevant  for 

z  <  0,  and  need  no  additional  consideration. 

The  technique  used  for  evaluating  the  kz  integration  consists  in  deform¬ 
ing  the  initial  contour  C  to  the  new  contour  C'  shown  in  Fig.  3.  Contour  C' 
begins  at  Re  kz  +  and  is  forced  to  pass  through  the  saddle  point  at  i2z/d|j2. 
The  contour  passes  below  the  kc  pole  (pole  to  the  left  of  contour)  and  finally 
goes  to  Re  kz  +  "  by  passing  above  the  ks  pole  (pole  to  the  right  of  contour). 
The  first  step  in  evaluating  the  kz  integral  (defined  by  I)  consists  of 
separating  the  individual  contributions  of  each  pole  by  writing  the  integrand 
in  terms  of  partial  fractions,  i.e.. 


(kz2-ks2rl(kz2-kc2)-l  ,  <ks2-kc2)-'  [Ckz2-ks2,->  -  <kz2-kc2)-l]. 


The  desired  integral  is  then  1=  Is  +  Ic,  where  Is  is  the  integral  having  the 
pole  at  kg  and  Ic  is  the  integral  having  the  pole  kc. 

The  Is  integral  can  be  readily  evaluated  because  the  pole  at  ks  remains 
far  from  the  saddle  point  for  all  values  of  kx.  Consequently,  the  influence 


of  the  pole  on  the  evaluation  of  the  saddle  point  contribution  Is  negliglible 
for  this  case.  Using  Cauchy's  theorem  and  approximating  the  path  contribution 
along  C'  by  the  standard  saddle  point  value  Isp  results  in 

Is  -  ISp  +  2tt1  Residue  (k2  =  ks).  (17) 

The  evaluation  of  Ic  using  the  C'  contour  requires  proper  attention  to 
the  possibility  that  the  pole  at  kc  may  be  found  close  to  the  saddle  point. 
However,  we  may  write  the  result  symbolically  in  the  form 


Ic  =  I'gp  +  2iri  Residue  (k2  *  kc),  (18) 

where  now  I'Sp  represents  the  saddle  point  result  properly  corrected  for  the 
presence  of  the  pole.  The  details  leading  to  the  evaluation  of  Ic  are  given 
in  Appendix  A;  it  suffices  to  mention  here  that  the  result  involves  the 
familiar  plasma  dispersion  function  evaluated  at  a  parameter  C  *  kcd|/2  - 
iz/d| . 

Combining  the  details  involved  in  evaluating  the  expressions  in  Eqs.  (17) 
and  (18)  transforms  Eq.  (15)  into  the  following  well  behaved  integral: 


-A  .  .  2 


Ep.^P.z)  =  r  djdx 


/ 


dki  kI2  J 1  (kj.p) 


[ki4(l-e1/e|1)2  +  4k04  (l-ki2/k02e  | )  ]  lf2 

exp  (-kj.2dx2/4)  [Gp,^  +  Kp>(^  +  QPj(j,]  ,  (19) 


gP,4>  =  ^  exP  (-ks2dB  2/4  +  iksz)  hp^  (kz  *  kg) 


(20) 


=  -(1/2)  exp  (-z2/d|2)  hp>(j,  (kz  =  kc)  Z( 4 ) 


(21) 


-12- 


Qp*<|)  =  ~exp  (-z2/dj2)  jdp>(j)  (kcd  g  +  12z/dn)/d|2 

+  ( 1  /2)hp, (j,(kz  =  i2z/dj2)  [ (i2z/d j )(z2/d | 2  +  ks2d|2/4)-1 

+(kcd  h  /2  +  lz/d,)-M}  (22) 

where  6p  =  1,  6^  =  0,  4  =  kcdn/2  -  iz/d || ,  and 

fm  -t2 

Z(4)=  /7  /  dt  +  10  ^  e'c2  .  (23) 


a  = 


0  Im  O  0 

1  Im  5  =  0 

2  Im  C  <  0 


(24) 


is  the  plasma  dispersion  function.  When  0  =  l(for  Im  5  =  0)  in  Eq.(24),  the 
integral  in  Eq.(23)  is  to  be  interpreted  as  Cauchy's  principal  value  integral. 

Physically,  the  terms  contained  in  Gp>(j,  represent  the  fields  associated 
with  the  shear  Alfven  mode.  For  u>/u)cj  «  1  this  terra  is  the  dominant  contrib¬ 
ution  to  the  radiation  pattern.  The  radial  dependence  of  the  shear  Alfven  mode 
is  primarily  determined  by  the  factor  kj^2Jj (k^p)exp(-k^2dx2/4)  in  the  inte¬ 
grand;  it  shows  that  the  shear  mode  fields  vanish  at  p  *  0  and  peak  at  some 
fraction  of  dj.,  i.e.,  they  create  a  ring-like  pattern  highly  collimated  along 
the  flux  tube  passing  through  the  hot  spot.  The  Kp>(^  terra  contains  the  con¬ 
tribution  of  the  compressional  mode,  which  for  u>  <  o>cj  is  much  smaller  than 
the  shear  mode.  The  fact  that  a  portion  of  the  spectrum  with  kj_  >  k^  does  not 
contribute  to  a  propagating  mode  (Fig.  2)  is  properly  taken  care  of  by  the  Z 
function  in  the  Kp  ^  term.  The  weak  excitation  of  the  compressional  mode  is 
evident  from  the  coupling  coefficient  hp  of  Eq  (15).  Finally,  the  Qp>4>  term 
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term  contains  purely  near  field  effects  associated  with  the  local  conduc¬ 
tivity  modulation. 

Ill .  Radiation  Pattern 

To  obtain  detailed  quantitative  information  about  the  spatial  variation 
of  the  electric  fields  produced  by  the  conductivity  modulation,  the  kj^  Inte¬ 
gration  in  Eq.  (19)  is  performed  numerically  for  specific  parameter  values. 

Figure  A  shows  the  radial  profile  of  the  magnitude  of  each  of  the 
electric  field  components,  scaled  to  |a|  -  ( A°z/0zo)mEo,  for  different  axial 
positions.  The  parameters  used  are  c/v^  *  390,  u>/u>cj  *  0.23,  Vj/u>  »  2xl0~4, 
k^d f  ■  k^dj^  *  1.  Figure  Aa  shows  the  radial  field  Ep  which  Is  dominated  at 
large  z  by  the  contributions  due  to  the  shear  mode.  The  characteristic  head¬ 
light  pattern  with  null  field  on  axis,  common  to  other  field  aligned  excita¬ 
tion  scenarios,3  exhibits  a  peak  at  k^P  ~  0.7  and  decays  rapidly  for  k^p  >  2 
(i.e.,  it  behaves  as  exp  (-p2/d^2)).  Figure  Ab  displays  the  azimuthal  field 
E^  which  is  much  smaller  than  Ep  (note  factor  of  10  in  scale).  In  addition  to 
the  headlight  pattern  associated  with  the  shear  mode  contribution,  E^  exhibits 
a  small,  but  slowly  decaying,  feature  at  large  p  corresponding  to  the  radiated 
compre8sional  mode. 

The  axial  dependence  of  Ep  and  E^  Is  shown  in  Fig.  3  for  a  fixed  radial 
position  near  the  peak  of  the  headlight  pattern,  i.e.,  k^p  *  0.7.  It  is 
evident  from  the  presentation  that  within  k^z  ■  1.5  the  shear  mode  attains  its 
asymptotic  value  with  a  polarization  I E p |  »  |E^|. 

The  individual  contributions  to  the  total  Ep  (i.e.,  arising  from  the  Gp, 
Kp,  and  Qp  terras  in  Eq.  (19))  are  shown  in  Fig.  6  as  a  function  of  radius  for 
different  axial  positions.  In  Fig.  6a  one  observes  that  the  shear  mode 
contribution  is  essentially  constant  beyond  k^z  «  0.5.  This  behavior  is  to  be 
contrasted  with  Fig.  5,  where  it  is  seen  that  the  total  Ep  does  not  reach  its 
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asymptotic  value  until  k^z  *  1.5.  The  difference  between  these  two  figures  is 
attributable  to  the  contributions  arising  from  the  near  field,  shown  In  Fig. 
be.  For  k^z  <  l  the  near  field  contribution  tends  to  cancel  the  pure  shear 
mode  contribution  and  results  in  the  net  pattern  shown  in  Figs.  4  and  5.  The 
radially  spreading  nature  of  the  compress ional  mode  is  visible  in  Fig.  6b  where 
it  is  clear  that  at  large  axial  values  k^z  =2.5  the  peak  moves  further  away 
from  the  magnetic  flux  tube  formed  by  the  shear  mode  of  Fig.  6a. 

The  previously  discussed  field  patterns  pertain  to  a  spherically 
symmetric  hot  spot,  i.e.,  dj  =  d^.  The  effect  of  changing  dj/dj^  on  Ep  is 
shown  in  Fig.  7  for  a  cigar-like  spot,  i.e.,  longer  along  the  field  line  than 
across  it  (k^d^  =0.5,  k^dj  =  1),  a  spherical  spot  (k^dj_  *  1,  k^dj  =  O,  and  a 
pancake-like  spot,  i.e.,  larger  perpendicular  to  the  field  than  along  it 
(k^d^  =  l,  k^dn  =  0.5).  Fig.  7  illustrates  that  the  headlight  pattern  due  to 
the  shear  mode  becomes  narrower  and  of  larger  intensity  for  the  cigar-like 
spot  than  for  the  spherical  one.  For  the  pancake-like  spot  the  headlight 
pattern  is  similar  to  that  obtained  for  the  spherical  spot,  but  the  magnitude 
of  the  electric  field  is  reduced. 

Another  quantity  of  interest  is  the  time  averaged  Poynting  vector  = 
(c/8n)  Re  (J3  x  _B*),  which  in  the  present  problem  has  a  natural  scaling  factor 
SQ  =  (c2/v^)  |a|2/8h.  The  radial  dependence  of  the  various  components  of  is 
shown  in  Fig.  8  for  selected  values  of  z,  and  for  the  same  conditions  leading 
to  Figs.  4-6.  It  is  clear  from  Figs.  8a  and  8b  that  at  these  low  frequencies 
power  is  radiated  predominantly  along  the  confining  magnetic  field  (z  direc¬ 
tion)  with  a  highly  collimated  headlight  pattern  passing  through  the  hot  spot. 

A  much  smaller  fraction  of  the  total  radiated  power  is  carried  away  by  the 
compressional  mode  out  in  the  direction  perpendicular  to  the  magnetic  field. 

In  addition,  there  is  a  circulating  Poynting  vector  component  in  the  azimuthal 
direction  shown  in  Fig.  8c.  This  circulating  component  does  not  lead  to 


net  radiation 


Finally,  it  is  useful  to  Introduce  a  scaled  figure  of  merit,  that  charac¬ 
terises  the  radiation  efficiency  of  this  system.  Of  the  possible  scaling 
factors  that  can  be  considered,  perhaps  the  more  meaningful  is  the  amount  of 
power  going  into  the  ohmic  modulation  at  frequency  to,  i.e.,  the  volume 
integral  of  A oz  EQ2/2.  Using  this  definition  leads  to 

// S  .  da 

n  *  ~ P - =2  ,  (25) 

IIIj  AczE0dv 


where  n  is  a  scaled  geometric  factor  that  must  be  evaluated  for  a  particular 
shape  of  the  hot  spot.  The  order  of  magnitude  of  n  is  Sz/S.  The  important 
parameter  dependences,  however,  are  contained  in  the  factor  multiplying  n  in 
Eq.  (26).  The  enhancement  due  to  v^  shows  up  in  the  form  (c/v^)2,  and  the 
expected  low  inefficiency  inherent  to  all  low  frequency  schemes  is  measured  by 
(u)/ozo).  Since  the  process  considered  is  linear,  the  dependence  on 
(Aoz/ozo)m,  i.e.,  the  depth  of  modulation,  enters  linearly.  Equation  (26)  may 
also  be  expressed  as  n  *  (ui/u)cj)( ve/ujce)(Aoz  /°zo)n,  which  is  another  way  of 
looking  at  the  parameter  dependence  of  the  efficiency,  and  showing  that  it  is 
density  Independent. 

IV.  Concl uslons 

The  radiation  pattern  generated  by  a  field  aligned  current  due  to  a 
localized  modulation  of  the  electrical  conductivity  at  a  frequency  u>  below  the 

ion  cyclotron  frequency  o)cj  has  been  calculated.  The  pattern  consists  of  a 
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shear  Alfven  mode,  a  compresslonal  Alfven  mode,  and  near  fields.  The  shear 
mode  propagates  in  a  highly  collimated  pattern  along  the  magnetic  field  lines 
passing  through  the  hot  spot  giving  rise  to  the  modulation,  while  the 
compresslonal  mode  is  emitted  over  a  broad  range  of  angles  and  moves  radially 
away  from  the  hot  spot.  For  u)/<uci  <  1,  nearly  90%  of  the  radiated  power 
appears  in  the  form  of  the  shear  mode  with  the  remaining  10%  being  taken  away 
by  the  compressional  mode. 

The  radiation  efficiency  for  this  system  is  found  to  scale  according  to 
(c/vA)2  (w/ozo)  (Aoz/oz0),  thus  indicating  the  enhancement  associated  with  the 
high  index  of  refraction  found  in  a  magnetized  plasma.  However,  the  inherent 
inefficiency  associated  with  the  generation  of  low-frequency  electric  fields 
in  a  highly  conducting  plasma  mitigates  the  enhancement  by  the  factor  a>/ozo. 

For  parameters  characteristic  of  the  auroral  ionosphere,  where  naturally 
occurtng  field  aligned  currents  are  present,  *  it  is  found  that  n  ~  10-5  and 
the  magnitude  of  the  transverse  electric  fields  associated  with  the  shear  mode 
can  be  comparable  to  the  magnitude  E0  of  the  field  driving  the  current  system. 
Although  this  low  frequency  generation  scheme  appears  inefficient  for  practi¬ 
cal  usage,  because  the  shear  mode  is  highly  collimated,  it  may  be  possible  for 
properly  placed  detectors  to  sample  this  component.5 

It  is  worth  noting  that  although  we  have  been  primarily  motivated  to 
consider  the  modulation  of  field  aligned  currents  arising  due  to  an  external 
electron  heating  source,  the  ideas  presented  here  and  the  techniques  for  calcu¬ 
lating  the  fields,  described  in  Sec.  II,  can  be  applied  to  situations  in  which 
field  aligned  modulations  arise  by  other  means.  For  instance,  in  a  recent 
experiment  in  a  research  tokamak  Borg,  et.al.6  have  modulated  the  current  in  a 
field  aligned  external  wire.  Preliminary  measurements  have  shown  the  genera¬ 
tion  of  a  shear  mode  having  a  headlight  pattern  similar  to  the  type  discussed 
in  this  paper.  It  is  also  found  by  these  experimentalists  that  the  radiated 


compressional  mode  is  much  weaker  than  the  shear  mode.  Another  situation  In 
which  related  phenomena  occur,  albeit  without  direct  external  control.  Is  In 
the  appearance  of  the  well  known  sawtooth  oscillations  attributed  to 

'l 

temperature  relaxations  In  tokamaks. 
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The  Integral  of  Eq.  (18) 

Ic  =  I'Sp  +  2tt1  Residue  (kz  =  kc)  (Al) 

contains  a  saddle  point  Integration  I’gp  where  the  analyticlty  condition 

Imposed  on  the  integrand  Is  affected  by  the  presence  of  a  first  order  pole  at 

kz  =  kc  c^at  is  in  the  neighborhood  of  the  saddle  point  at  kz  =  ksp. 

Following  the  treatment  of  Banos,2  this  situation  Is  properly  handled  by 

expressing  the  saddle  point  Integration  as 

I'sp  “  Xa  +  V  (A2) 

where  we  have  separated  the  regular  analytic  part  of  the  integral  Ia  from  the 
part  containing  the  effect  of  the  pole  Ip. 

From  Eqs .  (13)— (18),  we  write  the  integrand  of  Eq.  (Al)  as 

k~  g(k  ) 

®(kz)  “  v  2  J  f  2  exP  (— kz2d| 2/^  +  ikzz),  (A3) 

Kz  Kc 

where  g(kz)  is  an  analytic  function  In  the  complex  kz  plane.  Using  partial 
fractions  Eq.  (A3)  is  rewritten  as 

*(kz)  -  y  (k  .  +  g(kz>  exp  (-kz2d12/4  +  ikzz)  -  *_  +  *+.  (A4) 

The  first  term  of  Eq.  (A4),  <&_,  has  a  pole  at  kz  »  -kc  away  from  the  saddle 
point  (for  positive  z  (Fig.  3)).  Therefore,  its  contribution  to  the  saddle 
point  integration  is  calculated  in  the  usual  manner.  The  second  term  ♦+  with 

the  pole  at  kz  -  kc  passing  in  the  vicinity  of  the  saddle  point  is  properly 


-19- 


handled  by  expanding  the  analytic  function  g(kz)  about  kc,  which  yields 


MM-  1  +  «'(M  + 


z/  2  1  kz-kc 


+....] 


exp  (-kz2dj2/4  +  ikzz). 


where  the  primes  denote  derivatives  with  respect  to  kz. 

From  Eqs .  (A2),  (A4),  and  (A5)  we  can  write  the  Integrand  of  Ia  as 

1  g(k  )  g"(k  ) 

$a^kz^  =  J  l-k  +k  "  +  S'^kc^+  — JT~~  (kz-kc>  +  •••]  exP  (“kz2d|2/4  +  lkz2>»  (A&) 
z  c 

and  the  integrand  of  Ip  as 

1  s(k  ) 

*p<kz)  *  2  k'~k  exp  (-kz2dl2/4  +  lkz2>*  (A?) 

z  c 

The  task  of  evaluating  Ic  along  the  contour  C’  of  Fig.  3  is  now  greatly 
facilitated  by  rewriting  Eq.(Al)  as 


Ic  ■  Ia  +  lb  =  *a  +  [*p  +  2wl  Residue  (kz  *  kc>] 


From  Eq.  (A6)  it  Is  clear  that  the  integral  Ia  Is  easily  evaluated  along  the 
contour  C'  using  the  common  saddle  point  method,  while  the  terra  in  square 
brackets  in  Eq.  (A8)  containing  the  effect  of  the  pole  is  written  as 


lb  =  Ip  +  2ni  Residue  (kz  ■  kc) 


g(kc)  exp  (-z2/d|2)  Z(0, 


where  C  *  (kcd|/2  -  iz/d|)  and  7,(0  is  the  plasma  dispersion  functi< 
defined  in  Eqs.  (23)  and  (24). 


When  kc,  which  is  a  function  of  k^  (Fig.  3),  is  away  from  ksp,  the  value 
of  Z(C)  is  such  that  the  pole  has  negligible  contribution  to  the  saddle  point 
integration.  As  kc  gets  closer  to  ksp,  Z(c)  takes  values  that  properly 
account  for  the  presence  of  the  pole  near  kgp.  This  result  is  valid  for  all 
values  of  C  in  the  complex  kz  plane. 
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Fig.  I 


Figure  Captions 


Schematic  representation  of  the  model  showing  the  collimated 
headlight  pattern  of  the  radiated  shear  Alfven  wave. 


Fig.  2  Dispersion  relation  showing  kz2  versus  k^2  for  fixed  u)/<ocj  =  0.23 
and  nA  =  390. 

Fig.  3  Location  of  the  poles  and  saddle  points  in  the  complex  kz  plane  for 
the  kz  integral  of  Eqs.  (17)  and  (18)  (for  z  >  0).  The  initial 
contour  C  is  deformed  to  C'. 

Fig.  4  Radial  profile  of  the  normalized  magnitude  of  the  (a)  radial  and 
(b)  azimuthal  electric  field  components  at  two  axial  positions 
kAz=  0.5  (near  field)  and  kAz  =  2.5  (far  field),  w/o»cl  =  0.23, 
nA  =  390,  kAd#  =  kAdx  =1. 

Fig.  5  Axial  profile  of  the  normalized  magnitude  of  the  electric  field 
components  Ep  and  E^  at  a  fixed  radial  position  kAP  =  0.7.  Same 
parameters  as  in  Fig.  4.  Constant  magnitude  at  large  axial 
distance  is  due  to  the  shear  Alfven  mode. 

Fig.  6  Radial  profiles  of  the  components  of  Ep  at  the  same  axial  locations 
shown  in  Fig.  4  for  the  same  parameters:  (a)  shear  Alfven  wave; 

(b)  compressional  Alfven  wave;  (c)  near  fields. 

Fig.  7  Dependence  of  Ep  on  the  source  size  for  w/u>c£  =  0.23  and  nA  =  390. 

Heated  volume  is  (a)  cigar-like;  (b)  symmetric;  (c)  pancake-like. 

Fig.  8  Normalized  Poynting  flux  components  for  the  fields  of  Fig.  4. 

Over  90%  of  the  power  is  radiated  in  the  z  direction. 
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